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ABSTRACT 


A  stochastic  analysis  of  the  Supply  Corps  officer 
manpower  system  is  presented.  The  analysis  is  based  on  the 
concept  of  diffusion  approximations.  The  state  of  a  rank 
is  represented  as  the  superposition  of  a  deterministic  process, 
and  a  random  noise  variable,  both  of  which  are  derived 
utilizing  a  set  of  formulated  recursion  equations.  Numeri¬ 
cal  results  are  presented  for  a  3-rank  system  and  compared 
against  a  simulation  of  a  Markovian  manpower  model. 

A  Diffusion  Approximation  analysis  is  performed  utilizing 
United  States  Navy  Supply  Corps  data,  in  an  effort  to 
determine  the  probability  of  maintaining  each  rank,  any 
two  ranks,  and/or  the  entire  system  between  specified  limits. 
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I .  INTRODUCTION 


This  paper  attempts  to  examine  manpower  problems  uti¬ 
lizing  stochastic  techniques.  The  point  of  departure  is  the 
deterministic  theory  of  control  given  in  Bartholomew  [Ref.  1], 
and  the  object  is  to  see  how  concepts  and  strategies  origi¬ 
nating  in  that  context  perform  in  a  stochastic  environment 
of  the  kind  likely  to  be  encountered  in  practice.  Conclusions 
reached,  both  in  this  thesis  and  in  Bartholomew  [Ref.  3], 
suggest  that,  although  useful,  deterministic  analysis  is 
incomplete  since  it  fails,  for  example,  to  provide  required 
information  concerning  the  variance  of  individual  ranks. 
Without  such  information,  manpower  planners  may  attempt  to 
"fine  tune"  the  system  by  increasing  or  decreasing  promotion 
or  recruitment  rates  when  in  reality  the  fluctuations  are 
due  to  the  uncertainty  in  the  system. 

This  study  concentrates  its  efforts  on  one  particular 
manpower  system — that  being  the  Supply  Corps  of  the  United 
States  Navy.  An  effort  is  made  to  develop  a  stochastic  model 
that  provides  the  capability  for  simple  computation.  In 
addition,  it  will  give  the  mean  number  of  individuals  in 
each  rank,  as  do  deterministic  models;  and  as  added  features, 
will  provide  second  order  moments  and  distribution  of  force 
information.  Armed  with  these  key  elements  we  can  then 
attempt  to  answer  questions  related  to  the  maintenance  of 
desired  force  structures  and  with  the  ability  to  attain  such 
structures . 
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The  numbers  of  people  in  particular  ranks  change  over 
time  as  a  result  of  attrition,  promotion  and  recruitment. 


Some  of  these  flows  are  subject  to  management  control  while 
others  vary  due  to  a  multitude  of  circumstances  and  may  best 
be  modelled  stochastically.  Factors  such  as  the  need  to  offer 
adequate  career  prospects,  while  at  the  same  time,  maintain¬ 
ing  prescribed  force  levels  are  the  manpower  planners  chief 
goals.  A  problem  arises  in  that  the  planner  must  determine 
whether  these  goals  can  be  achieved  simultaneously,  and,  if 
so,  how  best  to  achieve  them. 

The  manpower  control  problem  has  received  a  good  deal  of 
attention  in  recent  years — see,  for  example,  Bartholomew 
[Refs.  1,2,3],  Davies  [Refs.  4,5]  and  Grinold  and  Stanford 
[Ref.  6] . 

Although  most  of  the  work  on  manpower  control  has  been 
based  on  stochastic  models,  the  analysis  has  been  determinis¬ 
tic  in  the  sense  that  only  expected  values  for  attrition  and 
promotion  rates  are  used.  More  realistically,  attrition  and 
other  flows  do  vary  randomly  and  are  therefore  better  repre¬ 
sented  by  random  variables.  Recent  work  which  does  attempt 
to  incorporate  randomness  into  manpower  control  are  Bartholomew 
[Refs.  2,3]  and  Kim  [Ref.  7]. 

In  this  thesis  we  will  assume  that  only  promotion  rates 
and  recruitment  rates  are  subject  to  control.  While  it  may 
be  true  that  attrition  can  be  dramatically  increased,  uti¬ 
lizing  techniques  such  as  "Reductions  in  Force”,  it  is  the 
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problem  of  being  unable  to  hold  down  attrition  rates  that 
plagues  the  present  day  Navy.  As  a  result,  the  model  developed 
in  this  thesis  utilizes,  as  an  input  parameter,  the  attrition 
rate  and  is  capable  of  determining,  in  a  deterministic  sense, 
the  promotion  and  recruitment  rates  required  to  maintain  the 
system  at  some  specified  level.  Henceforth,  the  concept  of 
being  able  to  maintain  a  specified  force  structure  will  be 
referred  to  as  maintainability;  i.e.,  a  structure  is  main¬ 
tainable  if  and  only  if  there  are  promotion  and  recruitment 
rates  such  that  for  each  rank  in  the  system,  the  expected 
number  of  individuals  who  depart  during  each  time  period  is 
equal  to  the  expected  number  of  individuals  who  arrive  during 
that  same  time  period.  The  majority  of  this  thesis  will  deal 
with  the  maintainability  of  desired  force  structures,  but  it 
turns  out  that  this  topic  cannot  be  treated  adequately  with¬ 
out  considering  the  question  of  whether  those  structures  can 
be  attained  from  present  levels. 

Since  for  all  practical  purposes,  recruitment  into  the 
Supply  Corps  takes  place  only  into  the  lowest  rank,  promotion 
control  alone  is  left  as  a  means  of  controlling  the  force 
distribution.  Control  by  promotion  has  the  disadvantage  that 
fluctuations  in  the  promotion  rate  are  liable  to  create 
internal  stresses  within  the  organization.  On  the  other 
hand,  it  was  shown  in  Bartholomew  [Ref.  1,  Chapter  4],  that 
in  a  deterministic  environment,  promotion  control  is  capable 
of  maintaining  a  much  wider  range  of  structures  than  recruit¬ 
ment  control  alone. 
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Having  discussed  the  various  areas  of  consideration,  the 
general  scope  of  this  paper  will  be  to  present  the  following: 


A.  A  basic  deterministic  Markovian  model  will  be 
derived  in  Chapter  II,  to  be  used  as  a  benchmark  with  which 
to  compare  the  stochastic  version. 

B.  Realizing  that  a  deterministic  representation  is  not 
representative  of  the  random  factors  that  affect  most  pro¬ 
cesses,  a  more  realistic  stochastic  model  will  be  derived 

in  Chapters  III  and  IV  utilizing  a  Diffusion  Approximation. 

C.  Results  of  a  simulation  of  a  stochastic  Markovian 
manpower  model  will  be  presented  in  Chapter  V  and  compared 
to  results  of  the  Diffusion  Approximation  as  developed  in 
this  thesis.  Shortfalls  of  the  approximation  will  be  dis¬ 
cussed  and  advantages  of  the  approximation  pointed  out. 

D.  A  comparison  of  the  deterministic  and  stochastic 
models  will  be  made  in  Chapter  VI.  In  particular,  the  corre¬ 
lation  between  adjacent  ranks,  and  its  effect  on  the  ability 
to  maintain  desired  force  structures  will  be  discussed. 

E.  Areas  of  possible  further  study  will  be  briefly 
mentioned  in  Chapter  VI.  An  effort  is  made  to  identify 
possible  research  topics,  whose  results,  when  used  in  conjunc¬ 
tion  with  the  model  derived  in  this  thesis,  may  reflect  to 

an  even  greater  degree  the  stochastic  nature  of  a  manpower 
system. 


TO.  ,  V - r 
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II.  DETERMINISTIC  ANALYSIS 


A.  A  DETERMINISTIC  MANPOWER  MODEL  FOR  THE  SUPPLY  CORPS 

Before  one  attempts  to  model  any  system,  one  must  be 
familiar  with  its  everyday  operations.  In  an  effort  to 
compare  the  results  of  the  approximation  against  a  benchmark, 
we  will  use  our  familiarity  of  the  system  to  make  several 
educated  assumptions  as  to  the  values  of  promotion  and  attri¬ 
tion  rates  within  the  supply  Corps.  A  compilation  of  these 
educated  guesses,  will  define  a  Markovian  transition  matrix; 
the  elements  of  which  are  the  fractions  of  individuals  in 
each  rank  that  stay  in  the  rank  or  are  promoted  in  each  time 
period.  This  matrix  will  then  be  compared  against  one  where 
historical  data  is  used  as  a  basis  for  determining  the  ele¬ 
ments  of  the  Markovian  transition  matrix.  Comparison  between 
the  two  will  answer  the  question  of  whether  organizational 
assumptions  made  are  reasonable. 

The  organizational  assumptions  that  we  will  make  are  as 
follows ; 

1.  We  will  assume  that  the  fraction  of  Supply  Corps 
officers  with  a  particular  number  of  years  in  a  specified 
rank  remains  almost  constant;  that  is,  we  will  assume  that 
the  fraction  of  the  supply  Corps  force  that  are  second  year 
Lieutenants  or  third  year  Lieutenant  Commanders,  etc., 
remains  the  same. 

2.  Under  present  promotion  patterns,  promotion  from 
Ensign  (q-1)  to  Lieutenant  Junior  Grade  (0-2)  occurs  at  the 
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2  year  point  in  an  officer's  career.  In  addition,  promotion, 
for  all  practical  purposes,  is  assured  for  all  Ensigns  com¬ 
pleting  their  second  year.  It  is  therefore  assumed  in  light 
of  assumption  1,  that  approximately  50%  of  all  0-1 's  are 
completing  their  second  year  and  will  therefore  be  promoted 
to  0-2. 

3.  Similarly,  it  takes  two  years  to  transition  between 
ranks  0-2  and  Lieutenant  Commander  (0-3);  however,  only  95% 
of  those  in  the  second  year  in  rank  0-2  get  promoted.  There¬ 
fore,  this  implies  that  50%  of  all  0-2' s  are  in  the  promotion 
zone,  and  95%  of  those  eligible  are  promoted.  This  informa¬ 
tion  yields  the  promotion  fraction  equal  to  0.475  (0.50  x  0.95). 

4.  Using  similar  logic  the  following  table  values  are 
computed : 

Years  Required  Assured  per-  Assured  Prenotion 
in  grade  prior  centage  of  promotion  fraction  of 


to  being  in 
zone 

total  grade 
in  zone 

rate  for 
in  zone 

all  members 
in  grade 

Ensign  (0-1) 

2 

0.5000 

1.0000 

0.5000 

Lieutenant 
Junior  Grade 
(0-2) 

2 

0.5000 

0.9500 

0.4750 

Lieutenant 

(0-3) 

5 

0.2000 

0.9000 

0.1800 

Lieutenant 

Conmander 

(0-4) 

7 

0.1429 

0.8000 

0.1143 

Oomnander 

(0-5) 

5 

0.2000 

0.6000 

0.1200 
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In  order  to  complete  the  transition  matrix,  it  is  necessary 


to  determine  reasonable  attrition  rates.  Attrition  rate  data 
was  obtained  from  historical  data  files  at  the  Naval  Military 
Personnel  Command.  Attrition  rates  cannot  be  forecasted  with 
great  accuracy,  due  in  great  measure  to  their  correlation  to 
many  economic  and  social  factors.  Thus  it  appears  reasonable 
to  use  a  simple  moving  average  as  a  predictor.  We  will  use 

vi(t-i)  +  Vi(t) 

2 

as  a  forecasted  value  for  the  current  attrition  rate  for 
stage  i  at  time  t+1. 

This  average  appears  to  be  a  better  forecasting  technique 
than  that  of  simply  using  the  last  attrition  rate.  A  reason 
for  not  considering  more  terms  in  the  attrition  forecast, 
is  the  fact  that  we  want  the  forecasted  attrition  rate  to 
reflect  only  the  recent  past  and  not  to  consider  forces  which 
may  have  affected  the  manpower  system  at  one  time  but  are  no 
longer  present. 

The  remaining  element  required  in  the  transition  matrix 
is  the  continuation  rate.  Since  the  sum  of  the  continuation 
rate,  the  attrition  rate  and  the  promotion  rate  must  equal 
one,  it  is  a  simple  matter  to  compute  the  remaining  elements 
of  the  transition  matrix.  It  should  be  mentioned  once  again, 
that  these  assumptions  are  not  related  to  any  historical 
data,  except  for  attrition  rates  and  existing  "in  zone" 
promotion  rates.  The  resulting  matrix  is  computed  merely  to 
compare  against  the  matrix  derived  in  the  next  subsection. 
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B.  PROMOTION  RATES  REQUIRED  TO  MAINTAIN  THE  SUPPLY  CORPS 

AUTHORIZED  FORCE  SIZE 

The  objective  of  this  subsection,  is  to  find  out  if  the 
force  levels  which  the  supply  Corps  has  been  authorized, 
are  in  fact  maintainable.  The  authorized  force  levels  men¬ 
tioned  are  included  in  Appendix  A,  and  are  obtained  from  the 
Officer  Programmed  Authorizations,  Fiscal  Years  1979-1983 
[Ref.  9].  In  an  effort  to  work  with  a  single  force  level, 
the  average  force  level  was  computed  for  the  five  years  men¬ 
tioned  above;  and  it  is  this  level  that  we  attempt  to  maintain. 
It  should  be  noted  that  over  the  5  years  under  consideration, 
neither  the  individual  ranks  nor  total  force  strength  vary 
significantly,  so  the  average  force  appears  to  be  a  reasonable 
force  level  to  maintain. 

If  we  determine  that  the  desired  force  can  be  maintained, 
we  will  then  have  a  transition  matrix  which  can  be  compared 
against  the  matrix  P^  derived  in  subsection  II. A.  This  will 
give  us  an  idea  of  how  drastically  the  system  would  have  to 
change  in  order  to  maintain  the  desired  force  levels. 

Before  we  proceed,  a  few  definitions  for  the  above: 

Q  =  a  six  comporent  row  vector  whose  elements 
are  the  number  in  each  rank. 


P  =  a  6  x  6  matrix  whose  elements  contain  promo¬ 
tion  and  continuation  rates.  (See  matrix 
P,,  subsection  II. A  for  the  general  form 
or  this  matrix.) 

T 

V  =  a  column  vector  whose  elements  contain  the 
appropriate  attrition  rates  for  each  rank. 
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R  =  the  recruitment  vector,  which  directs  the 
flow  of  recruits  into  the  proper  ranks 
(fraction  of  recruits  for  each  rank). 

In  the  Supply  Corps,  the  recruitment  vector  is  not  gen¬ 
erally  amenable  to  control,  i.e.,  recruitment  for  the  most 
part  occurs  into  the  lowest  rank,  i.e.,  R=  [1,0, 0,0, 0,0]. 

In  such  cases  we  wish  to  know  what  degree  of  control  can  be 
exercised  by  the  promotion  rates.  For  a  fixed  size  organiza¬ 
tion,  what  we  wish  to  do  is  to  find  a  transition  matrix  P, 
which  allows  us  to  maintain  our  desired  structure  Q.  This 
implies  that  P  must  satisfy 

Q  -  QP  +  Q  VT  R  (II.B.l) 

T 

where  Q  V  R  yields  the  number  of  recruits  necessary  to 
exactly  equal  the  number  of  attritions  at  each  time  period. 

As  was  noted  in  subsection  II. A,  the  matrix  P  is  super¬ 
diagonal  in  form,  i.e.,  non-zero  elements  only  in  the  main 
diagonal  and  the  diagonal  above  it.  Bartholomew  [Ref.  1] 
shows  that  when  P  has  the  super-diagonal  form,  equation  (II .B.l) 
admits  a  unique  solution  which  may  be  tested  to  see  whether  its 
elements  are  non-negative.  If  the  matrix  P  is  found  to  be 
non-negative,  then  the  structure  Q  can  be  maintained. 

As  was  mentioned  previously,  while  we  assume  that 
recruitment  takes  place  only  into  the  lowest  rank,  we  recog¬ 
nize  the  fact  some  recruitment  does  occur  in  ranks  0-2  and 
0-3.  However,  the  number  recruited  into  these  ranks  is  small 
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and  hence  will  be  neglected.  Also,  the  attrition  rates 
are  as  considered  previously,  i.e., 

V  =  [0.0272,0.1355,0.1121,0.0390,0.0790,0.1622]. 

We  note  that  equation  (II.B.l)  is  equivalent  to 

QP  =  Q-QVTR 

and,  at  this  point  it  becomes  more  convenient  to  abandon 
the  matrix  notation  temporarily  and  to  proceed  with  the  com¬ 
ponents  as  does  Bartholomew  [Ref.  1]: 

°ipn  *  °i  -  Ri 

k 

°1P12  +  °2P22  31  °2  "  R2  ^Vi 

k 

°4P56  *  °4P66  *  °6  “  R6  J^^i  * 

(The  above  are  components  of  the  respective  vectors  or 
matrices  mentioned  at  the  beginning  of  this  subsection.) 

Eliminating  P..  using  P..  =  1  -  V.  -  P.  ...  (i  *  1,2,..., 

XI  11  jL  1/1*1 

k-l)  and  using  our  special  structure  for  the  vector  R  yields 
the  following  equation: 
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(II. B. 2) 


Pi, i+1 


k 

l  Q-:  V  ./Qi 
j=i+l  3  ]  1 


In  other  words,  equation  (II.B.2)  says  that  the  proportions 
required  to  be  promoted  from  grade  i  must  equal  the  number 
leaving  from  all  grades  above  grade  i  divided  by  the  size 
of  grade  i. 

For  example,  given  a  point  estimate  for  the  attrition 
rate,  and  given  that  the  promotion  rates  for  0-6' s  is  essen¬ 
tially  zero;  it  is  possible  to  compute  the  number  of  0-6* s 
who  depart  the  system.  Given  this  value,  one  can  then 
compute  the  number  of  0-5's'who  must  be  promoted  each  year 
to  exactly  maintain  the  required  level  of  0-6's.  The  process 
continues  in  this  fashion  until  we  compute  the  required  num¬ 
ber  of  new  accessions  required  to  maintain  the  0-1  Rank.  It 
is  in  this  way  that  promotion  rates  can  be  calculated  to 
exactly  maintain  the  desired  structure  (assuming  that  we  are 
at  the  structure  we  wish  to  maintain) .  If  we  are  at  some 
other  force  level  and  we  wish  to  attain  our  goal,  similar 
calculations  can  be  done,  but  review  must  be  done  to  ensure 
that  the  resulting  promotion  rates  are  within  acceptable 
bounds  as  specified  by  organization  policy. 

For  example,  suppose  we  wish  to  maintain  the  following 
structure: 
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Total  Force  Level  =  3955 


Rank 

Strength 

Attrition  Rates 

o 

i 

M 

621 

0.0272 

0-2 

625 

0.1355 

0-3 

995 

0.1121 

0-4 

863 

0.0390 

0-5 

617 

0.0790 

0-6 

234 

0.1622 

We  first  compute  the  number  of  0-6' s  lost.  This  is  equal  to 
(234)  (0.1622)  =  37.9548.  To  maintain  the  desired  structure 
37.9548  0-5 's  on  the  average  must  be  promoted.  Therefore, 
the  promotion  rate  required  is  (37.9548)  ~  (617)  «  0.0615. 

An  important  note  to  make,  is  that  promotion  rate  is  used 
in  relation  to  the  total  rank  strength,  to  determine  what 
the  rate  would  be  in  regards  to  "In  the  promotion  zone" 
strength,  one  must  first  determine  the  number  to  be  promoted, 
and  then  divide  that  number  by  the  number  "in  the  zone".  In 
our  example,  promotion  rate  related  to  "In  the  zone"  would  be 
37.9548  f  (number  of  0-5' s  in  zone).  The  other  promotion 
rates  cam  be  similarly  computed,  and  as  a  final  result  we 
obtain  the  matrix  P2. 

Comparing  this  matrix  with  reveals  that  the  reasonable 
guesses  for  the  existing  system  are  very  close,  in  most  cases, 
to  the  rates  necessary  for  maintaining  the  desired  system  on 
the  average.  This  implies  that  very  little  has  to  be  done  in 


TRANSITION  MATRIX 


way  of  changing  promotion  rates  if  we  simply  want  to  get 
to  and  maintain  the  desired  levels.  It  should  be  stressed, 
that  once  a  specified  force  level  has  been  shown  to  be 
maintainable,  it  can  be  attained  from  any  starting  force 
level  as  was  stated  in  Ref.  1.  In  the  Supply  Corps,  present 
rank  levels  are  not  equal  to  desired  rank  levels.  However, 
if  promotion  rates  were  modified  to  those  suggested  by  P2, 
we  would  eventually  reach  a  system  of  manpower  stocks  whose 
expected  value  in  each  rank  would  equal  the  desired  levels. 

A  technique  that  suggests  a  method  to  speed  up  the  process 
of  convergence  from  present  levels  to  desired  levels  will  be 
discussed  in  a  later  chapter. 
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III.  THE  DIFFUSION  MODEL 
The  system  we  wish  to  model  is  as  depicted  below: 


The  variables  of  interest  are  as  follows: 

Q  =  the  force  available  for  recruitment, 

°  assumed  to  be  large. 

a(t)  5  the  accession  rate  at  time  t 

u.  (t)  =  the  promotion  rate  for  all  members  in 

rank  i  at  time  t 

v.  (t)  =  the  attrition  rate  for  all  members  in 

rank  i  at  time  t 

Q. (t)  i  the  Random  Variable  which  is  the  force 
level  in  rank  i  at  time  t. 

Other  variables  of  interest  to  be  incorporated  into  the 
model  are: 

q.(t)  =  the  expected  fraction  of  the  total  force 

1  in  rank  i  at  time  t  (computed  as  fraction 

of  original  total  force  level) 


E[Qi(t)l 
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p 

W.  (t)  =  Promotion  random  noise  variable, 

1  distributed  N(0,1)  (when  i  *  0, 

W  p(t)  relates  to  accession) 
o 

W.L(t)  =  Loss  random  noise  variable,  distributed 

1  11(0,1) 

X^(t)  =  standardized  random  noise  variable. 

The  probability  that  a  member  of  rank  i  gets  promoted 
at  time  t  is  equal  to  u^(t);  this  implies  that  the  proba¬ 
bility  that  the  member  does  not  get  promoted  is  equal  to 
l-u^(t).  We  will  assume  that  whether  or  not  an  individual 
is  promoted  is  independent  of  others  in  the  system  at  time 
t.  Thus,  the  number  promoted  will  have  a  binomial  distribution 
with  expected  value  u^(t)Q^(t)  and  a  variance  of 
u±(t)  [1  -  ui(t) ]Qi(t) . 

Using  similar  reasoning  the  number  of  attritions  in  rank 
i  at  time  t  is  assumed  to  have  a  binomial  distribution  with 
an  expected  value  of  v^(t)Q^(t)  and  a  variance  of 
vi(t) [1  -  vi(t) ]Qi(t) . 

In  rank  1  we  have  a  special  case,  since  we  must  consider 
the  accession  rate  as  a  function  of  the  total  force  level. 

We  will  assume  that  there  is  a  large  population  of  potential 
recruits  of  size  Qq.  We  assume  the  number  of  people  that 
enter  the  system  at  time  t  has  a  binomial  distribution  with 
mean  a(t)QQ  and  a  variance  of  a{t)  [l-a(t)]QQ.  In  what 
follows  we  will  usually  take  Qq  to  be  the  initial  force 
strength. 
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If  one  deals  with  a  system  which  is  sufficiently  large, 
one  can  invoke  the  Central  Limit  Theorem.  Therefore  if  Qq 
is  "large",  this  would  imply  that  Q^(t)  is  "large"  and 
therefore  the  distribution  of  number  promoted  can  be  approxi¬ 
mated  by  a  Normal  Distribution  with  mean  u^(t)Q^(t)  and  a 
variance  of  u^(t)  [1  -  u^(t) ]Q^(t) .  Similar  results  apply  to 
the  number  of  attritions  in  each  rank. 

Now,  recall  if  X  is  distributed  normal  with  mean  y  and 
2  2 

variance  a  [X  -  N(u,<j  )].  Then 

aX  +  b  -  N(ap  +  b,a^a^) . 


Let 


WiP(t)  -  N ( 0 , 1)  , 


then 


Qi(t)ui(t)  +  ^Qi  (t)  ui  (t)  [1  -  u.(t)  ]  W^U) 

-  N(Qi(t)ui(t),Qi(t)ui(t) [1  -  ui(t)]) 

which  is  a  normal  distribution  with  mean  and  variance  equal 
to  the  mean  and  variance  specified  for  the  binomial  distribu¬ 
tion  of  the  number  promoted  at  time  t. 
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Similarly , 


Qi(t)vi(t) 


+  /Q~(t)vTTt)  [1  -  vi(t)  ]  W^tt) 


-  N(Qi  (t)  vi  (t)  ,Qi  (t)  vi  (t)  [1  -  vi  (t)  ]  ) 


is  a  normally  distributed  random  variable  with  mean  and 
variance  equal  to  the  mean  and  variance  specified  for  the 
binomial  distribution  of  the  number  attrited  at  time  t. 

For  each  interval  of  time  (t,t+l]  we  can  define  the 
change  in  rank  i,  i.e.,  Qi(t+1)  -  Q. (t),  as  a  combination 
of  normally  distributed  random  variables. 


Qi(t+1)  -  Qi(t)  * 

Ui-lQi-llt)  +  l-1  -  ]Qi_1(t)  W^/lt) 

(promotions  into  rank  Q^) 

-  ui(t)Qi(t)  -  ✓ui  (t)  [1  -  u^t)  iQ^t)  WiP(t) 

(promotions  out  of  rank 

-  vi(t)Qi(t)  -  ✓vj.tt)  [1  -  vi(t)  ]Qi(t)  WiL(t) 

(attritions  out  of  rank  Q^) 


From  this  relationship,  Q^t+1)  can  be  solved  for  in  terms 
of  Q^(t).  Hence,  we  assume  for  i  >  1, 


Qi(  t+1) 


/ui-l ^  "  ui-l ^Qi-1 (t) 


Wi-1 


(t) 


-  /u^t)  [1  -  ui(t)  ]Qi(t)  WiP(t) 

-  ^vi(t)  [1  -  v..(t) ]Qi(t)  WiL(t) 

+  [1  -  u±  (t)  -vi(t)]  Qi(t)  (III.l) 

where  (W^P(t),  i  =  0,1,..., r-1}  and  (W^L(t),  i  =  l,...,r} 
are  independent  random  variables,  having  a  normal  distribution 
with  mean  zero  and  variance  equal  to  one.  In  addition  all 
the  W  random  variables  are  independent  of  (Q^(t),  i  =  l,...,r}. 
Let 


Xi(t) 


Qi (t)  -  Qoqi(t) 


/Q- 

o 


(III. 2) 


If  Qq  is  sufficiently  large,  then  X(t)  =  [X^ ( t) , X2 (t) , . . . , Xr (t) ) 
should  have  approximately  a  multivariate  normal  distribution 
with  mean  zero. 

From  equation  (III. 2)  we  have 


Qi(t)  *  Qoqi(t)  +  Xi(t)  1 

In  Appendix  B  we  replace  the  variables  Qi(t)  in  equation 
(III.l)  by  the  expression  in  (III. 3)  and  then  take  the  limit 
as  Qq  -►  ».  This  procedure  yields  recursive  equations  for 
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{ qi  ( t: ) ,  i  =  l,...,r}  and  {Xi(t),i  =  The 

{qi(t),i  =  l,...,r}  are  constants  and  {X^(t),i  =  l,...,r} 

is  a  multivariate  normal  random  variable  with  mean  zero. 

Hence  in  light  of  (III. 3),  Q(t)  will  also  have  a  multivariate 

normal  distribution,  at  least  in  the  limit  as  Q  «.  In 

o 

the  next  chapter  we  will  present  a  recursive  scheme  to 
evaluate  qi(t),  Var[X^(t)]  and  cov [X^ (t) , X^ ( t) ] . 
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IV.  THE  RECURSIVE  SCHEME 


Only  results  will  be  presented  in  this  section.  Deriva¬ 
tions  of  the  results  will  be  presented  in  Appendix  B  of  this 
paper . 

The  following  are  definitions  required  for  the  recursions 
which  follow: 


Var[Xi(t)]  =  itk  ( t) 

Cov [X. (t) ,X . (t) 1  =  m.  .  (t) 

i  j  1  >  J 

RESULT  1:  A  RECURSIVE  SCHEME  FOR  COMPUTING  q^t),  t  =  1,2,... 
i  ■  1:  q1(t+l)  *  a(t)  +  [1  -  u1  (t)  -  v1  (t)  ]  q^t) 

i  >  1:  qi(t+l)  =  ui_1  (t)  qi_1  (t)  +  [1  -  u±  ( t)  -  (t)  ]  q^t) 

RESULT  2:  A  RECURSIVE  SCHEME  FOR  COMPUTING  VarlX^t)],  t  =  1,2,... 

i  *  1:  m1(t-tl)  =  [l-u1(t)  -v1(t)]2  m^(t) 

+  a (t) [1  -  a(t) ] 

+  (t)  [1  -  ux  (t)  ]  +  v1  (t)  [1  -  v1(t)  ]  Jq^t) 
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i  >  1:  m^t+l)  =  [u^ft)]2 

+  [1  -  ui(t)  -  vi(t) ] 2  mi(t)  +  ui_1 [1  -  ui_1(t) ]qi_1(t) 

+  {u± (t)  [1  -  ui(t) ]  +  Vi(t)  [1  -  vi(t) ] }qi(t) 

+  2u .  . (t)  [1  -  u. (t)  -  v. (t) ]m.  ,  . (t) 

1  — i.  X  X  X  X  t  X 

RESULT  3:  A  RECURSIVE  SCHEME  FOR  COMPUTING  COV [X. (t) , X . (t) ] , 
t  =  1,2,...  J 

i  /  1  and  j  =  i+1: 

m.  . (t+1)  *  u.  . (t)u. (t)m.  ,  • (t) 

X# j  X- X  X  X  X*X 

+  u.  ,  (t)  [1  -  u. (t)  -  v . (t) ]m.  ,  . (t)  + u. (t) 

i-i  j  j  i  i,J  A 

•  [1  -  vk  (t)  -  vi  (t)  (t) 

+  [1  -  u. (t)  -  v. (t) ] [1  -  u. (t)  -  v. (t) ]m.  .(t) 

i  i  j  j  j 

-  uA(t)  [1  -  u^t)  ]qi(t) 
i,j  ft  i  and  j  r1  i+1 

+  [1  -  Uj  (t)  -  Vj  (t)  j 

+  (t)  [1  -  \xi  (t)  -  vA  (t)  (t) 

+  11  -  ( t)  -  v^t)  ]  [1  -  Uj  (t)  -  Vj  <t)  ]mi;^  (t) 


— r 


Yr 
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i  =  1  and  i  /  2 


m1  j (t)  =  uj_x (t) [1  -  ux (t)  -  vx (t) 

+  [1  -  (t)  -  v1  (t)  In^ ^  (t) 

It  should  be  noted  that  since  the 

Cov[X1(t)  ,X2(t)  ]  =  E[X1(t)X2(t)  ]  -  E[X1(t)]E[X2(t)  ] 

and  since  X^(t)  is  distributed  N(0,a2);  the  Cov [X^ ( t) ,X2 (t) ] 

is  equal  to  the  mixed  second  moment,  designated  m^2(t). 

Similarly,  the  variance  of  Xi(t)  is  equal  to  the  second  moment 

of  Xi(t)  =  E[Xi2(t) ]  =  mi(t) . 

Given  q^(t)  and  m^t)  we  have  the  distribution  of  Qi(t), 

since  each  Q^(t)  is  distributed  Normally  with  a  mean  equal 

to  Q  a- (t)  and  a  variance  equal  to  Q  m. (t) .  As  a  result  of 
OX  ox 

the  possible  correlation  between  ranks,  of  greater  importance 
to  us  is  the  Multivariate  Normal  Distribution  of  Q(t).  No 
routine  was  available  that  provided  the  capability  of  com¬ 
puting  the  probability  of  maintaining  all  the  ranks  within 
certain  bounds;  however,  a  bivariate  normal  distribution  was 
used  to  compute  probabilities  of  maintaining  adjacent  ranks 
within  plus  and  minus  one  standard  deviation  from  the  expected 
rank  level.  The  reason  for  desiring  this  capability  will 
be  discussed  in  Chapter  VI. 
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V.  BINOMIAL  SIMULATION  VS.  THE  DIFFUSION  APPROXIMATION 


Prior  to  relying  on  any  approximation  of  a  real  world 
system,  a  user  should  have  a  feel  for  how  close  the  approxi¬ 
mation  comes  to  reflecting  the  true  system.  In  this  chapter 
we  will  simulate  the  year  to  year  operation  of  an  imaginary 
Markovian  3-rank  system.  The  system  will  be  operating  with 
promotion,  attrition  and  accession  rates  such  that  the  expected 
value  of  each  rank  is  constant,  i.e.,  the  system  is  main¬ 
tainable  at  some  specified  level.  The  starting  parameters 
are  listed  below,  both  in  tabular  and  graphical  form. 


Qi<t) 

#  of  individuals 

ui(t) 

Promotion  Rate 

v^t) 

Attrition  Rate 

Rank 

1 

1000 

0.70 

0.10 

Rank 

2 

2000 

0.10 

0.25 

Rank 

3 

1000 

0.00 

0.20 

The  accession  rate  is  specified  as  0.20  of  an  available  body 
of  recruits  estimated  at  4000  individuals  each  year. 

Graphically  the  system  would  appear  as  depicted  below: 
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The  numbers  in  the  parentheses  are  the  expected  number  of 
individuals  flowing  across  the  appropriate  arcs  during  each 
time  period.  Once  again  we  note,  that  in  a  deterministic 
sense,  the  system  is  in  exact  balance,  i.e.,  the  number 
leaving  each  rank  is  exactly  equal  to  the  number  flowing  into 
the  rank. 

Let 

Qq  i  total  number  in  the  system 

Qi(t)  =  number  of  individuals  in  rank  i  at  time  t 

P.(t)  i  number  of  individuals  promoted  from  rank 

i  at  time  t 

W. (t)  r  number  of  individuals  attrited  from  rank 
i  at  time  t 

u^(t)  =  promotion  rate  for  rank  i  at  time  t 

v^(t)  =  attrition  rate  for  rank  i  at  time  t 

a(t)  =  accession  rate  at  time  t. 

Note  that  P^t)  and  W^(t)  are  not  independent  random  variables 
since 

Pitt)  +  Wi(t)  <  Qi ( t) 

A  program  BIMOD,  outlined  in  Appendix  C,  was  designed  to 
simulate  this  system.  BIMOO  operates  as  follows: 


1.  The  program  is  designed  to  replicate  the  operation  of 
a  manpower  system,  for  a  user  specified  number  of  times. 
(Maximum  number  of  replications  allowed  is  1000.)  In  the 
results  which  will  follow,  the  number  of  independent  repli¬ 
cations  used  was  300.  A  replication  consists  of  35  years  of 
promotions,  accessions  and  attritions  of  the  manpower  system 
under  consideration.  The  time  increment  in  each  block, 
corresponds  to  "t"  in  all  the  variables  mentioned  above. 

The  reason  blocks  of  35  years  were  chosen  for  comparison 
with  the  approximation,  results  from  the  fact  that  the 
diffusion  approximation  reaches  a  constant  covariance  matrix 
at  the  35  year  point.  By  terminating  each  replication  at  35 
years,  and  storing  the  value  of  each  rank  at  the  end  of  each 
run,  we  can  then  compute  an  estimate  for  the  covariance  matrix 
of  the  simulation  at  the  35  year  point. 

2.  For  each  time  period  t,  BIMOD  generates  a  binomial 
random  deviate  with  parameters  Q^t)  and  [u^t)  +v^(t)]. 

This  deviate  represents  the  total  number  of  individuals 
leaving  rank  i  at  time  t;  both  through  promotion  and  through 
attrition.  Call  this  random  variable  L^(t),  i  =  1,  . ..,  r. 

3.  For  each  time  period,  we  also  generate  independent, 
binomial  random  deviates  with  parameters  L^(t)  and 
u^(t)/[U^(t)  +v^(t)J,  i  -  l,...,r;  and  Qq  and  a(t).  These 
deviates  represent  the  total  number  of  promotions  from  rank 
i  at  time  t,  P.(t);  and  the  total  number  of  accessions  into 
the  system  at  time  t,  A(t) . 
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4.  The  Q^(t)  are  updated  for  each  time  period  utilizing 
the  following  recursive  equations 

Q-L  (t+1)  =  Q1(t)  -  LL(t)  +  A(t) 

Qi(t+1)  =  Qi(t)  -  Li(t)  +  P^U)  for  1  =  2,3 

In  comparing  the  simulation  against  the  diffusion  approxi 
mation,  we  would  like  to  reduce  the  comparison  to  the  random 
terms.  We  do  this  in  order  to  determine  how  good  the  approxi 
mation  is  in  predicting  the  means,  the  variances,  and  the 
covariances  of  the  random  terms. 

As  was  discussed  in  Chapter  III,  equation  (III. 3)  is  the 
cornerstone  of  the  diffusion  approximation 

Qi{t)  =  Qoqi(t)  +  (IV. 1) 

We  could  describe  the  system  of  our  simulation  in  similar 
terms, 

Mi(t)  =  Qomi(t)  +  RA(t)  (IV. 2) 

where  QQm^(t)  is  the  expected  value  of  rank  i  at  time  t, 

Mj (t)  is  the  observed  value  of  rank  i  at  time  t,  and  R^(t) 
is  the  random  residual  term  for  rank  i  at  time  t.  The  com¬ 
parison  of  (IV. 1)  and  (IV. 2)  would  commence  with  an  examina¬ 
tion  of  the  R^(t)  residual  terms 


(IV. 3) 


R^t) 


M. (t)  -  Q  m. (t) 
1 _ O  1 


To  see  if  the  F^(t)  do  in  fact  exhibit  an  approximate  normal 
distribution  with  mean  equal  to  zero. 

BIMOD  was  run  with  the  following  input  values: 


NUMBER  OF 

REPLICATIONS 

300 

YEARS  CONSIDERED  IN  EACH  REPLICATION 

35 

INITIAL  FORCE  LEVELS 

Q1 

1000 

°2 

2000 

°3 

1000 

PROMOTION 

RATES 
^  (t) 

0.70 

u2  (t) 

0.10 

u3(t) 

0.00 

ACCESSION 

RATES 

a(t) 

0.20 

ATTRITION 

RATES 

vx(t) 

0.10 

v2(t) 

0.25 

v3(t) 

0.20 

At  the  end  of  each  35  year  replication,  the  value  of 
each  rank  was  stored,  and  the  residual  values  were  computed 
as  in  equation  (IV. 3).  The  expected  values  used  in  equation 
(IV. 3)  were  the  initial  force  levels  as  stipulated  above, 
since  they  are  the  levels  we  expect  to  maintain.  When  the 
300th  replication  was  completed,  the  mean  value  of  the  residuals 
in  each  rank,  R^,  was  computed,  and  this  in  turn  was  used  to 
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compute  the  estimate  for  the  covariance  matrix  of  the  35th 
year  residuals.  The  formula  used  in  computing  the  covariance 
matrix  is 


Cov(Ri , R j ) 


=  l 


(Ri-Ri)(Ri-V 

n-2 


where  the  R^  are  the 

stored  values 

at 

the  end  of  each  repli 

cation. 

The 

following  is 

a  compilation 

of 

the  results  obtained 

from  the 

simulation. 

E[Ri(t) ) 

Kurtosis 

Skewness 

Rank  1 

0.0398 

2.7468 

0.2054 

Rank  2 

0.0162 

2.6670 

0.0248 

Rank  3 

0.0041 

2.8470 

0.0708 

Kurtosis  and  skewness  information  was  obtained  from  a 
histogram  of  the  total  run  residuals.  It  should  be  noted 
that  the  Kurtosis  Coefficient  for  a  normally  distributed 
random  variable  is  equal  to  3.0  and  the  skewness  should 
equal  0. 

This  first  level  of  comparison  of  the  R^t)  terms  with 
the  X^(t)  terms  seems  to  indicate  that  the  assumptions  of 
normality  with  mean  equal  to  zero  does  in  fact  provide  a 
reasonable  approximation,  with  a  possible  exception  of  rank 
1,  as  indicated  by  its  skewness  coefficient. 
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The  next  level  of  comparison  between  R^(t)  and  X^(t) 
involves  the  predicted  variances  and  covariances.  Utilizing 
BIMOD  the  following  covariance  matrix  was  obtained  for  the 
R^(t)  residual  terms. 


Rank  1 

Rank  2 

Rank  3 

Rank  1 

0.2035 

0.0000 

-0.0022 

Rank  2 

0.0000 

0.4747 

0.0245 

Rank  3 

-0.0022 

0.0245 

0.2783 

and  this  is  compared  to  the 

following  X^(t) 

covariance  l 

obtained  from  the 

diffusion 

approximation  program  MNPWR 

outlined  in  Appendix  C. 

Rank  1 

Rank  2 

Rank  3 

Rank  1 

0.2448 

-0.0210 

-0.0005 

Rank  2 

-0.0210 

0.5059 

-0.0289 

Rank  3 

-0.0005 

-0.0289 

0.2373 

At  first  glance,  it  would  appear  that  the  approximations 
for  the  covariance  terms  did  fairly  well.  Of  particular  note, 
however,  is  the  fact  that  the  approximation  specifies  a 
negative  correlation  between  ranks  2  and  3,  while  the  simu¬ 
lation  specifies  a  positive  correlation  for  the  same  ranks. 

In  order  to  better  understand  what  is  going  on,  the 
variance  of  the  covariance  terms,  i.e.,  the  variance  of  the 

^R1 2  ^R1  ~Rj \  for  all  end  of  run  R. ,  was  computed  to 

n-2  1 
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determine  whether  the  diffusion  approximation  estimate  falls 
in  the  range  Cov(R^,Rj)  ±  Std  Dev [Cov(R^,Rj ) ] .  The  following 
are  the  variance  results. 


Ranks 

Cov  -  1  Std  Dev 

Simulation 

Covariance 

Cov  +  1  Std  Dev 

1,2 

-0.0173 

0.0000 

0.0173 

1,3 

-0.0162 

-0.0022 

0.0118 

2,3 

0.0047 

0.0245 

0.0443 

While  the  covariance  range  specified  above  still  implies 
a  positive  covariance  for  ranks  2  and  3,  the  covariance 
values  under  consideration  are  small  enough  to  imply  that 
the  approximation  estimate  for  the  covariance  is  quite  good. 

The  next  important  matter  of  note,  is  the  fact  that  in 
both  matrices  above,  the  covariance  terms  are  essentially 
zero;  and  the  fact  that  ranks  appear;  to  be  independent  of 
one  another  is  borne  out  by  very  low  correlation  coefficients 
between  ranks,  both  in  the  simulated  system  and  in  the  diffu¬ 
sion  approximation. 

A  possible  explanation  for  the  resulting  independence 
could  be  that  since  we  do  not  vary  promotion,  attrition,  or 
accession  rates  in  either  the  simulation  or  the  diffusion 
approximation,  the  probability  laws  of  the  models  are  the 
same  as  those  in  a  model  in  which  the  motions  of  individuals 
through  the  system  are  independent.  It  is  known  that  if 
the  input  random  variables  are  independent  and  have  Poisson 
distributions  and  the  motions  of  the  individual  in  the  system 
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are  independent,  then  the  number  of  individuals  in  stage  i 
at  time  t,  Y^(t),  are  asymptotically  independent  with  marginal 
Poisson  distributions  as  t  +  *  [Refs.  1,8]. 

In  comparing  the  variances,  we  note  that  the  diffusion 
approximation  provides  reasonable  estimates  for  the  variance 
of  the  simulation  residuals. 

In  conclusion,  one  could  possibly  simulate  a  6-rank  sys¬ 
tem  and  obtain  better  results  for  the  variances  and  covariances 
of  the  ranks.  However,  due  to  the  enormous  number  of  deviate 
generations  required  to  produce  representative  output,  it 
is  felt  that  the  much  simpler  computations  of  the  diffusion 
approximation  far  outweigh  their  loss  of  accuracy. 


38 


VI.  DIFFUSION  APPROXIMATION  RESULTS 


In  this  final  chapter,  we  will  present  the  results 
obtained  when  the  6-rank  Supply  Corps  system  is  modeled 
utilizing  the  Diffusion  Approximation  Program,  listed  in 
Appendix  C.  In  addition  a  program  called  MINI  will  be  dis¬ 
cussed  and  its  results  reintroduced  into  a  second  run  of 
the  Diffusion  Approximation  to  the  Supply  Corps  system. 

The  final  subsection  of  this  chapter  will  provide  some 
recommendations  for  further  study. 

A.  A  DIFFUSION  APPROXIMATION  APPLIED  TO  THE  PRESENT  LEVELS 
OF  THE  SUPPLY  CORPS 

In  Chapter  II,  we  determined,  at  least  in  the  determinis 
tic  sense,  that  under  specific  assumptions,  desired  rank 
levels,  as  defined  by  Ref.  9,  could  be  maintained.  We  now 
wish  to  utilize  those  same  attrition  rates 


V(t)  =  [0.0272,0.1355,0.1121,0.0390,0.0790,0.1622] 


in  the  Diffusion  Approximation  and  to  apply  them  to  present 
Supply  Corps  force  levels.  These  levels  are  presented  below 


Rank 

Level 

0-1 

632 

0-2 

691 

0-3 

961 

0-4 

792 

0-5 

616 

0-6 

177 
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These  values  were  obtained  from  the  Supply  Corps  Officer 
Personnel  section  of  the  Naval  Supply  Systems  Command,  and 
were  current  values  for  310X  officers  as  of  September  1979. 

The  program  MNPWR  listed  in  Appendix  C,  incorporates  the 
recursive  equations  derived  in  Chapter  IV.  A  feature  of 
this  program  is  that,  given  the  present  levels,  desired 
levels,  and  estimates  for  attrition  rates  specified  pre¬ 
viously,  the  program  will  compute  the  following  fixed  promo¬ 
tion  rates, 

U  =  [0.5099,0.3711,0.1210,0.1005,0.0615,0.0000] 

and  a  fixed  accession  rate  "a"  of  0.0862.  Note  that  these 
promotion  and  accession  rates  are  those  which  are  necessary 
to  attain  and  maintain  the  desired  force  levels  given  that 
we  start  at  the  present  force  levels,  and  are  the  same  as 
those  specified  in  Section  II. B  of  this  thesis. 

Therefore,  in  the  first  of  a  series  of  computations, 
we  attempt  to  see  how  long  it  will  take  to  attain  the  desired 
force  levels  using  the  promotion  and  accession  rates  specified 
above.  Note,  we  are  not  worried  about  maintainability,  since 
Ref.  1  states  for  the  deterministic  case,  that  if  a  structure 
is  maintainable,  it  is  also  attainable.  However,  in  the 
stochastic  sense  maintainability  has  a  slightly  different 
connotation.  In  considering  methods  of  control,  we  are  now 
interested  in  methods  that  not  only  provide  the  expected 
rank  levels,  but  also  provide  manpower  planners  estimates  for 


the  variance  of  each  rank.  Returning  to  the  initial  compu¬ 
tation,  the  reason  we  are  not  concerned  about  maintainability 
is  that  if  the  program  is  able  to  compute  non-negative 
accession  and  promotion  rates,  then  the  structure  is  main¬ 
tainable. 

A  program  run  was  made  to  determine  how  long  it  would 
take  to  attain  the  desired  force  structure.  At  a  point  50 
years  into  the  future,  the  expected  number  in  each  rank  was 
approximately  at  the  desired  level.  This  length  of  time  is, 
without  much  thought,  unacceptable.  A  technique  which  allows 
faster  convergence  to  the  desired  levels  will  be  described 
in  the  next  subsection. 

The  following  is  a  listing  of  the  expected  values  and  the 
covariance  matrix  of  the  system  at  time  50. 

EXPECTED  VALUE  OF  EACH  RANK  AT  TIME  50 


Rank 

E [Rank] 

Desired 

0-1 

620.99927 

621 

0-2 

624.99902 

625 

0-3 

994.99780 

995 

0-4 

862.98999 

863 

0-5 

616.91138 

617 

0-6 

233.85883 

234 
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Analysis  of  the  covariance  matrix  will  be  reserved  until 
the  next  subsection. 


B.  MINIMUM  CONVERGENCE  TIME 

The  quickest  method  of  attaining  the  desired  force  from 
the  present  force  would  be  simply  to  promote  (or  attrit) 
all  shortages  (or  excesses)  in  one  time  period.  In  the  case 
being  examined,  this  would  be  totally  infeasible  both  for 
legislative  and  for  practical  purposes.  To  attain  the 
desired  levels  in  one  time  period  would  involve  a  very  large 
number  of  officer  promotions,  which  would  be  totally  unaccepta¬ 
ble  from  the  organizational  standpoint. 

The  approach  taken  in  program  MINI  ( see  Appendix  C) ,  is  a 
Dynamic  Programming  Approach.  Given  what  are  reasonable 
bounds  for  promotion  and  accession  rates,  the  program  will 
minimize  the  difference  from  present  levels  to  desired  levels, 
at  each  time  period,  for  specified  attrition  rates. 

Utilizing  the  following  bounds  for  accession  and  promotion 
rates. 


Lower  Bounds 
(accession  and  promotion 
rates  required  to  maintain 
the  system) 

Upper  Bounds 
(110%  of  the  lower 
bound  rates) 

acces¬ 

sion 

0.0862 

0.0948 

0-1 

0.5099 

0.5609 

0-2 

0.3711 

0.4082 

0-3 

0.1210 

0.1331 

0-4 

0.1005 

0.1106 

0-5 

0.0615 

0.0677 

0-6 

0.0000 

0.0000 
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and  the  attrition  rates  specified  at  the  beginning  of  this 
chapter,  we  could  reduce  the  convergence  time  required  for 
the  expected  number  of  individuals  in  each  rank  to  attain 
the  desired  level  from  about  50  years  to  at  most  8  years. 

CONVERGENCE  TIME  REQUIRED 


Rank 

Time  Required 
for  Convergence 

0-1 

4  years 

0-2 

3  years 

0-3 

1  year 

0-4 

8  years 

0-5 

6  years 

0-6 

8  years 

As  an  output  of  program  MINI,  we  obtain  the  promotion  rates 
required  to  minimize  the  shortfall  at  each  time  period.  The 
following  is  a  listing  of  the  promotion  rates  suggested  by 
MINI. 

MINI  PROMOTION 
AND  ACCESSION  RATES 


TIME 

0-1 

0-2 

0-3 

0-4 

0-5 

ACCESSION  RATE 

1 

0.5099 

0.3902 

0.1331 

0.1106 

0.0677 

0.0862 

2 

0.5099 

0. 3754 

0.1331 

0.1106 

0.0677 

0.0862 

3 

0.5104 

0.3829 

0.1331 

0.1106 

0.0677 

0.0862 

4 

0.5286 

0.3904 

0.1331 

0.1106 

0.0677 

0.0891 

5 

0.5292 

0.3904 

0.1331 

0.1106 

0.0677 

0.0893 

6 

0.5292 

0.3904 

0.1331 

0.1088 

0.0677 

0.0893 

7 

0.5292 

0.3904 

0.1331 

0.1066 

0.0677 

0.0893 

8 

0.5179 

0.3790 

0.1260 

0.1015 

0.0622 

0.0875 

Note  that  between  years  7  and  8  the  promotion  and 
accession  rates  begin  to  decrease.  This  is  a  direct  result 
of  the  fact  that  we  are  approaching  the  desired  levels,  and 
following  year  8,  we  then  use  u  and  "a"  as  specified  in 
Section  VI. A. 

Utilizing  these  computed  accession  and  promotion  rates 
as  input  parameters  to  the  Diffusion  Approximation  yields 
information  on  how  the  system  would  react  to  system  controls. 
The  following  is  a  summary  of  the  Diffusion  Approximation 
results  at  the  8  year  point. 


Rank 

Desired 

Level 

Expected  Value 
Approximation  Result 

0-1 

621 

621.00903 

0-2 

625 

624.95386 

0-3 

995 

994.87744 

0-4 

863 

863.02441 

0-5 

617 

617.00854 

0-6 

234 

233.97964 

Of  even  greater  significance  are  the  variance/ covariance 
results.  Below  is  the  covariance  matrix  obtained  at  the  8 
year  point.  Variances  for  each  rank  are  on  the  main  diagonal. 

A  striking  result  becomes  evident  if  one  compares  the 
expected  value  and  the  variance  of  each  rank.  The  indepen- 
den  c  motions  property  is  still  evident,  even  with  promotion 
control,  i.e.,  the  expected  value  in  each  rank  is  approximately 
equal  to  the  variance  of  the  rank  (as  one  would  expect  for 
Poisson  distributed  random  variables) .  What  this  may  imply. 
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is  that  in  this  case  even  though  we  vary  the  promotion  and 
accession  rates,  the  changes  allowed  in  those  rates  are 
much  too  small  to  counteract  the  independent  motions  property. 
As  a  result,  bivariate  probabilities  taken  indicated  that 
the  ranks  are  essentially  independent. 

If  in  fact  the  ranks  are  relatively  independent,  one 
could  then  compute  expected  value  plus  and  minus  one  standard 
deviation  of  force  level  in  each  rank,  and  know  that  with  the 
assumption  of  a  normal  distribution  for  force  levels,  we 
have  approximately  a  68%  chance  of  being  in  that  interval. 

From  our  calculations  we  get  the  following  results  for  upper 
and  lower  bounds  of  force  level . 


Rank 

Lower 

Bound 

Desired 

Level 

Upper 

Bound 

0-1 

596 

621 

645 

<N 

1 

o 

598 

625 

651 

0-3 

962 

995 

1027 

0-4 

833 

863 

892 

0-5 

591 

617 

642 

0-6 

218 

234 

249 

Given  these  values,  manpower  planners  need  not  overreact 
to  random  deviations  in  force  levels  within  the  above  bounds. 
If,  for  example,  the  0-4  rank  were  to  fall  to  as  low  as  833, 
this  would  not  necessarily  imply  that  we  should  increase 
next  year's  accessions;  it  could  be  rationalized  as  a  reason¬ 
able  value  as  compared  to  the  variance  computed  by  the  Diffu¬ 
sion  Approximation.  Without  the  variance  estimates,  planners 


do  not  possess  the  tools  necessary  to  determine  whether  or 
not  fluctuations  in  each  rank  can  be  explained  by  the  internal 
randomness  of  the  System,  or  whether  the  fluctuations  are 
due  to  external  forces  which  require  changes  in  accession 
and/or  promotion  rates. 

C.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  Diffusion  Approximation  hinges  on  reasonable  estimates 
for  the  attrition  rate  of  each  rank.  Therefore,  Recommenda¬ 
tion  1  would  be  to  derive  a  useful  estimating  relationship 
for  the  forecasting  of  attrition  rates. 

Another  possibility  for  an  extension  of  this  thesis, 
would  be  to  consider  confidence  intervals  for  attrition 
rates,  and  to  input  these  intervals  with  the  Diffusion  Approxi¬ 
mation.  Analysis  of  the  results  obtained  would  be  interesting 
in  regards  to  the  independent  motions  observed  when  point 
estimates  were  used  for  attrition  rate. 

It  is  clear  from  the  foregoing  analysis  that  maintaining 
a  grade  structure  in  a  stochastic  environment  is  much  less 
straightforward  than  a  deterministic  analysis  of  the  problem 
would  suggest.  Under  rather  special  conditions  we  have  seen 
that  it  is  possible  to  maintain  a  structure  in  spite  of  random 
variation  in  each  rank.  With  better  forecasting  of  attrition 
rates,  it  is  felt  that  the  Diffusion  Approximation  will  prove 
to  be  an  invaluable  tool  to  manpower  planners. 
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APPENDIX  B 


DERIVATIONS 

B.l  DERIVATION  OF  THE  q.(t+l)  RECURSION,  THE  VARIANCE 
RECURSION,  AND  THE  RELATIONSHIPS  FOR  X^t) 

Q-^t+l)  -  Q^t)  =  a(t)Qo  +  /a(t)  [1  -  a (t) ]Q0WqP (t) 

-  u^tJQ^t)  -  (t)  [1  -  (t)  ]  Q1  (t)'  (t) 

-  V-^tjQ^t)  -  Jv1  (t)  [1  -  V1  (t)  ]Q1  (t)'  WxL(t) 

(for  definition  of  the  mentioned  variables,  see  Chapter  IV). 
Substituting  /Q^  X^ (t)  +  QQq1 (t)  for  Q^(t)  and  X^Ct+l) 

+  QQq1(t+l)  for  Q1(t+l),  we  have 

✓Q^lXj^t+l)  -  X1(t)  +  Vq^  qi(t+l)  -  qL(t)  ]  ■ 

a  (t)  +  /a(t)  [1  -  a  (t )  WQP(t) 

-  ux(t)  [X^t)  +  /Q^  qx(t)  ] 

-  V^tt)  [l-u1(t)][(l//5^)X1(t)  +  q1(t)l  WXP  (t) 

-  v1(t)  [X^Ct)  +  ^5”  q^lt)] 
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X1(t+1)  -  Xx(t)  = 

/Q^  qx(t)  -  ^  q1(t+l)  +  /Q^  a(t)  +  /a (t)  [1  -  a (t) ]  WQP(t) 

-  u-^t)  X^t)  -  u^t)  /Q^  q^t.) 

-  ^ux(t)  [1  -  u1(t)][Q./«/Q^)X1(t)  +qi(t)]  W^t) 

-  v1(t)  Xx(t)  +  v1(t)  /Q^  qj_(t) 

-  a/v  l(fc)  11  ‘  vl(t)  3  Kl/>/Q^)X1(t)  +  qx(t)  7  WxL(t) 

Heuristically ,  since  X^(t)  is  a  normal  random  variable, 

X1(t+1)  -x1(t)  must  be  finite  and  therefore  any  terms  with 

/Q^  must  approach  zero  as  Qq  •*  ®.  Therefore  in  the  limit  as 

Q  -*■  00 
O 


X.^t+1)  -  xx(t)  -  /a(t)  [1  -  a(t>  J  WoP(t) 

-  ux(t)  Xx(t)  -  /u^Tt)  [1  -  u]7tTJq^Tt)  W1P(t) 

-  vx(t)  Xx(t)  -  v^Tt)  [1  -  v1(t)]q1(t)  WxL(t) 


Rewriting,  we  have 


X1(t+1) 


[1  -  U1 (t)  -  (t) ]  X^ (t) 


+  /a(t)  [1  -  a(t)  ]  WQp(t) 

-  Ai”Tt)  [1  -  u1(t)  ] q1  ( t )  wxp(t) 

-  (t)  [1  -  (t)  ]q1  (t)  WxL(t)  (B.1.1) 

It  also  follows  that  in  the  limit  as  Qq  -  * 

q^  (t+1)  -  q^t)  =  a  (t)  -  (u^t)  +  v^tHq^t) 

or 

q.^t+1)  ■  a(t)  +  [1  -  u^t)  -  v^(t)  Jq^t)  (B.1.2) 

Using  a  similar  approach  for  Q^(t+1)  -  Q^(t),  i  f  1,  the 
following  result  was  obtained  in  the  limit  as  Qq  -*  ». 

Xi(t+1)  =  ui_1(t)Xi_1(t)  +  [1  -  ui(t)  -  vi(t)  ]  X^t) 

+  ATi_1(t)  [1  -  W^^tt) 

-  /ui(t)  [1  ~  uA(t)  Iq^t)  WiP(t) 

-  /vi(t)  tl  -  v.(t)]qi(tT  W.L(t)  (B.1.3) 
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and 


qi(t+l)  -  qi(t) 


ui-l  (t)  ~  Iu±(t)  +  vi(t)]qi(t) 


or 


qi(t+l)  =  ui_i (t) qi_1 ( t)  +  [1  -  ui(t)  -  vi(t) ]qi(t) 

(B.1.4) 

We  now  have  expressions  for  all  the  X^(t+1)  and  q^(t+l) 
in  terms  of  X^(t)  and  q^(t)  respectively. 

In  order  to  calculate  the  variances  and  covariances 
of  the  {X^ (t) , i =  1, . . . ,r } ,  we  must  first  determine  the  first 
and  second  moments  of  the  X^(t+1). 

EjX^t+1)]  =  [1  -  uL(t)  -  vx(t)  ]  E[XL(t)J 

+  /a(t)  [1  -  a ( t) ]  E[WoP(t)] 

-  ^u1 (t)  [1  -  u^ (t) ] q^ (t)  E[W1P(t)] 

-  /vx(t)  [1  -  vL(t) ]qx(t)  E [WXL (t) ]  (B.1.5) 

Recall  that  {W^P  (t)  ,i  =  0,1, . . .  ,r-l}  and  {W^L (t)  ,  i  =  1 , .  .  . , r } 
are  normally  distributed  with  a  mean  of  zero.  The  random 
variable  X^(t)  represents  the  randomness  in  the  system  at 
time  t.  Initially,  since  the  system  is  starting  from  a  known 


rm. . 
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set  of  rank  levels,  X^{0)  is  defined  to  be  zero.  As  a 
result,  equation  (B.1.5)  implies  that  E[X^(t+l)]  =  0  for  all 
time  t.  Further, 

Var[X1(t+l)]  =  [1  -  ux (t)  -  v1(t) ]2  Var [Xx(t) ] 

+  a(t)[l-a(t)]  Var  [WQP  (t)  ] 

+  u^t)  [1  -  u^tt)  lq^t)  Var[W1P(t)] 

+  (t)  [1  -  v1(t)  J q^ ( t>  Var[W1L(t)] 

Since  the  variances  of  the  VT  random  variables  are  equal 
to  one, 

VarfX^t+l)]  =  [1  -  ux  (t)  -  vx  (t)  ]  2  Var  [Xx(t)  ] 

+  a(t)[l-a(t)]  +  u^  (t)  [1  —  u^(t)  ]q^(t) 

+  v1(t)  [1  -  (t) ]q1(t) 

Since  the  expected  value  of  X1 (t)  is  equal  to  zero,  the 
variance  of  X^(t)  is  equal  to  its  second  moment  which  we 
will  call  m^(t).  Therefore  we  have  the  following: 

ra^t+l)  *  [1  -  ux(t)  -  vx(t)]2  m1(t)  +  a(t)  [1  -  a(t)  ] 

+  (t)  [1  -  u1  (t)  ]  +  (t)  [1  -  (t)  ]  >  q^  ( t) 
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Using  similar  reasoning,  we  get  a  general  recursive  scheme 
for  i  >  1: 


nu  (t+1 )  =  [ui_1(t)  ]2 

+  [l-ui(t)  -  vi(t)  ]2  m^t) 

+  (t)  [i  -  ui-;L  (t)  ]  q^ft) 

+  {u±(t)  [1  -  u± (t) ]  +  vi(t)  [1  -  vi(t) ] >q±(t) 
+  2{u.  .  (t) [1  -  u. (t)  -  v. (t) ] }  m.  ,  . (t) 

1—  X  1  X  1 

Note,  however,  for  i  >  1  we  get  a  covariance  term, 
the  next  section,  we  will  obtain  a  recursive  scheme  for 
computing  the  covariances. 

B.2  DERIVATION  OF  COVARIANCE  FORMULAS 

Recall : 


X^t+1) 


=  [1  -  ^  (t)  -  Vj.  (t)  ]  X1(t) 

+  /a(t)  [1  -  a (t) ]  WQP(t) 

-  (t)  [1  -  \i^  (t)  ]  q1  (t)  WxP(t) 

-  /v^t)  [1  -  v^t)  ] (t)  W^ft) 


In 


If  we  let 


b^t)  =  1  -  ui(t)  -  v^(t) 

c  (t)  =  a(t)  [1  -  a  (t)  ] 

(t)  =  ui(t) [1  -  ui(t) ] 
ki(t)  =  v^tHl  -  vi<t)J 

then 

X^t+1)  =  bL(t)  X1(t)  +  /cTt T  WQP(t) 

-  /h^Tt)q^TtT  W1P(t)  -  /k^TtT 

Similarly  recall 

X2(t+1)  =  ux(t)  Xx(t)  +  b2(t)  X2(t) 

+  /h1 (t) q1 (t)  W1P(t)  -  /h2 ( t) 

-  /ki (t) q2 (t)  W2L(t) 
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It  is  assumed  that  at  time  zero  the  Cov [X^ ( t) , (t) ]  =  0 
for  t  =  0, 

CtovtX^t+1)  fXjtt+l)  ]  ■  E[X1(t+l)X2(t+l)]  -  EtXj^tt+DJEt  Xj  (t+1)] 
x1(t+i)x2(t+i)  = 

[b-^tJX^t)  +/c(tr  WQP(t)  -  /h^tTq1  (t)  WLP(t)  -  /^(tjq^t)  (t)  ] 

•  tuL(t)X1(t)  +b2(t)X2{t)  +  ^1(t)q1(t)  WxP(t)  -  4i2(t)q2(t)  W2P(t) 

-  (t)q2(t)  W2L(t)] 

=  b1(t)u1(t)x12(t)  +b1(t)b2(t)xL(t)x2(t)  +b1(t)^h1(t)q1(t)  x1(t)W1P(t) 

-  b^(t)  /hj (t) q2 (t) Xj^  (t)  W2P (t)  -  bx(t)  /Tc2(t)q2(t)  X^tJW^t) 

+  uL(t)*/cTtT  x^tjw^t)  +  b2(t)  /STtT  x2(t)wQp(t) 

+  /cttjh^tjq^t)  W^tlW^t)  -  /c(t)h2(t)q2(t)  W2P(t)WQP(t) 

-  /C(t)k2(t)q2(t)  W2L(t)WQP(t) 

-  ^(tj^qitr^W  x^tjw^t)  -  b2(t)^h1(t)q1(t)  X^tW^t) 

-  h1(t)q1(t)  rwxP  (t)  ] 2  +  ^  (t) qx  (t) hj  ( t) q2  (t) W2P (t) (t) 

+  /hx  (t)  qx  (t)  k2  (t)  q2  (t)  W**  (t)  V*  (t)  -  u^t) /^(t)  (t)  Xx  (t)  W^t) 

-  b2(t)/k1(t)q1(t)  X^tW^t)  -  qx(t)  ^(td'k^t)  W1P(t)W1L(t) 

(t) q1'(t)h2(t)q2(t)  W2P  (t) WXL (t)  +  ^(tjq^tjk^tiq^t)  W^ttjW^t) 
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EtX^t+DX^t+l)]  = 

b±  ( t)  (t)  E  [X^2  (t)  ]  +  b1  (t)  b2  (t)  E  [3Cj_  (t )  ( t)  ]  +  ^  ( t)  /h]_(t)q1(t)  E  [3^  (t)  W.^  (t)  ] 

-  bL  ( t)  /h2  (t)  q2  (t)  E  [Xj.  (t) W2P  (t)  ]  -b1(t)/k2(t)q2(t)E[X1(t)W2L(t)] 

+  u^tJ^cTtTElX^tW^t)]  +  b2(t)v/cTtTE[X2(t)WoP(t)] 

Wc(t)h1(t)q1(t)  E [W1P  (t) WQP (t)  ]  -/c (t) h2  ( t) q2  (t)  E[W2P(t)WQP(t)  ] 

~^c  (t)  (t)  q2  (t)  E[W2L(t)WoP(t)]  -  ^(tJ^^Ttjq^Ttf  it)  ] 

-b2(t)  ^Ttj^tTE^ttJw/tt)!  -  h^tJq^tJEKW^t))2] 

+  /h1(t)q1(t)h2(t)q2(t) 

-  (t)  ^  ( t)  qx  (t)  E  [3^  ( t)  WXL  ( t)  ]  -  b2  ( t)  ^  ( t)  q1  (t)  E  [3^  ( t)  (t)  ] 

-  qx(t)  ✓h1(t)k:L  (t)  EtW^ttJW^tt)]  +/k1(t)q1(t)h2(t)q2(t)  E[W2P  (tJW^ft)  ] 

Since  the  X^(t)  are  assumed  independent  of  the  W^(t) 

P  2 

and  E[(W^  { t) )  J  =  1,  the  above  simplifies  to 

m12(t+l)  =  EtX^t+DX^t+l)  ]  *  b^t)  u^t)  m^t) 

+  bj^t)  b2(t)  n>12(t)  -  h1<t)  q^t) 
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to: 


m12*t+D  =  Cl  —  (t:)  —  V1  (t)  1  u^t)  m1  ( t ) 

+  [1  -  (t)  -  v1  (t)  ]  [1  -  u2  (t)  -  v2  (t)  ]  m12(t) 
-  ^  (t)  [1  -  ux  (t)  ]  q^  ( t) 

A  similar  approach  can  be  taken  for  all  cross  product 
terms,  and  the  following  general  relationships  will  be 
obtained: 

i  =  1,  j  >  2 


ml,j(t+i)  *  u.^tMl-u^t)  -v±{t))  mlfj_1 


(t) 


+  [1  -u,  (t)  -v,  (t)  ]  [1  -  u.  (t)  -  v.  (t)  ]  m,  .(t) 

1  1  J  J  ±9  J 


i  >  1,  j  =  i+1 


m.  . (t+1)  *  u.  . (t)  u. (t)  m.  ,  . (t) 

Xf  j 


+  u.  .  (t)  [1  -  u.  (t)  -  v.  (t)  ]  m  (t) 

i-i  j  j  i  i»J 

+  u^(t)  (1  -  uA(t)  -  v^t)  ]  mi(t) 

+  [1  -  u. (t)  -  v. (t) ] [1  -  u. (t)  -  v. (t) ]  m.  .  (t) 
13  33 


-  uA(t)  U  -  u^t)  ]  qi(t) 
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i>  j  >  1.  j  £  i+1 


m.  . (t+1)  =  u.  , (t)  u.  , (t)  m.  ,  .  ,  (t) 

X  #  J  IX  JX  IX, JX 


+  u.  . (t) [1  -  u. (t)  -  v. (t) ]  m.  ,  . (t) 

i-x  3  j 


+  u  .  (t)  (1  -  u.  (t)  -  v.  (t)  ]  m.  .  .  (t) 

J— X  X  1  x,j "X 


+  [1  -  ui  (t)  -  vi  (t)  ]  [1  -  Uj  (t)  -  Vj  (t)  ] 


To  fill  out  the  remainder  of  the  covariance  terms,  recall 


that 


nKjtt+l)  =  mj i  (t+D  • 


(t) 
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